Parametrically driven dark solitons 
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We show that unhke the bright sohtons, the parametricaUy driven kinks are immune from in- 
stabilities for all dampings and forcing amplitudes; they can also form stable bound states. In the 
undamped case, the two types of stable kinks and their complexes can travel with nonzero velocities. 
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The parametric driving is well known to be an effi- 
cient way of compensating dissipative losses of solitons 
in various media. Examples include surface solitons in 
vertically oscillating layers of water j^, ^; light pulses 
in optical fibers under phase-sensitive amplification |^ 
and in Kerr- type optical parametric oscillators [Q; mag- 
netisation solitons in easy-plane ferromagnets exposed to 
oscillatory magnetic fields in the easy plane Q . A serious 
problem associated with the parametric energy pumping, 
however, is that the driven solitons are prone to oscil- 
latory instabilities which set in as the driver's strength 
exceeds a certain — often rather low — threshold j|, || . 

With a few notable exceptions, the parametrically 
driven solitons considered so far had the form of pulses 
decaying to zero at spatial infinities. These were solu- 
tions of the nonlinear Schrodinger (NLS) equation with 
the "self-focusing" nonlinearity: 



iipt + ipxx + 2|V'PV' -ip = hip* - iji/j, 



(1) 



where ip is the amplitude of a nearly-harmonic station- 
ary wave oscillating with half the driving frequency, 7 
the damping coefficient, h the driving strength, and * 
indicates complex conjugation. However, in a number 
of applications the amplitude equation of the parametri- 
cally driven wave turns out to have the nonlinearity of 
the "defocusing" type: 



iipt + h->Pxx - IV'PV' + V' = hip* - ijTp. 



(2) 



In fiuid dynamics, the "defocusing" parametrically driven 
NLS (||) describes the amplitude of the oscillation of the 
water surface in a vibrated channel with a large width-to- 
depth ratio [|[ 0|. (On the contrary, the "focusing" equa- 
tion (|l|) pertains to the case of narrow channels.) The 
same equation (^) arises as an amplitude equation for the 
upper cutoff mode in the parametrically driven damped 
nonlinear lattices In the optical context, it was de- 
rived for the doubly resonant x''^'' optical parametric os- 
cillator in the limit of large second-harmonic detuning 
1^. Next, in the absence of damping, stationary solutions 
iP = My + iM, of eq.@ minimise the Ginzburg-Landau 
free energy for the anisotropic XY model, F = J J^dx, 
where 

= i(c>,M)2 - (1 + h)M^ + iM^ + 2hM^ + J^o, 

and M = {0, My, Mz). This model was used to study 
domain walls in easy-axis ferromagnets near the Curie 



point |]T^. Nonstationary magnetisation configurations 
were considered in the overdamped limit: ipt = —SF/6ip* 
pi] , [1^ . The damped hamiltonian dynamics ipt — 
—i6F/dip*—j'ip provides a sensible alternative; this is pre- 
cisely our eq.(|^). Finally, for 7 = there is yet another, 
independent, magnetic interpretation of eq.(|^); this will 
be introduced below. 

Localised structures characteristic of a defocusing 
medium are domain walls, or kinks, also known as "dark 
solitons" in the context of nonlinear optics. The purpose 
of this note is to explore the stability and bifurcations of 
the parametrically driven kinks and their bound states. 

Two stationary kink solutions of (H) are available in 
literature. One is usually called the Neel, or Ising, wall: 



''Pn{x) = iAia,'Dh{Ax)e 



Here ^ 1 + y/h'^ - 7^ and 6 



iarcsin^ 



(3) 
SI- 

For 7 = 0, the Neel wall coexists with the Bloch wall: 

Tpsix) = -iAta,-ah{2Vhx) ± Vl - 3ft.sech(2\//ix), (4) 

= I + h . Originally, the "magnetic" terminology 
was motivated merely by the fact that \tp\ = in the core 
of the wall (||), and so the point a; = is a phase defect 



similar to Neel points in solid state physics |12 . On the 



contrary, in the core of the Bloch wall the phase changes 
smoothly — and this is analogous to Bloch walls in fer- 
romagnets. Below we show that the analogy with mag- 
netism is in fact much deeper than originally thought. 
Letting h — ^ — 0, the Neel wall becomes the usual, un- 
driven, dark soliton whereas the Bloch wall degenerates 
into a flat solution. 

Both Bloch and Neel walls admit a clear interpretation 
in other physical contexts as well. For example, when 
eq.ra) is used to model the Faraday resonance in water 
[H, or chains of coupled pendula |H , both solutions de- 
scribe transitions between two domains oscillating 180° 
out of phase. The phase of the oscillation is discontinu- 
ous across the Neel wall and the amplitude changes over 
a narrow region; hence the wall appears as a highly lo- 
calised defect. Conversely, the Bloch wall has a smooth 
helicoidal structure, with the amplitude varying over a 
wider interval. It might therefore be tempting to ex- 
pect that in the region of their coexistence {h < |), the 
Neel wall should be unstable against the decay into the 
"smoother" (Bloch) wall, and this is indeed the case in 
the Ginzburg-Landau and Klein-Gordon counterparts of 
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eq.(|) (TT[ |T||. Surprisingly, the NLS dynamics turn 
out to be very different. 

We will show, numerically, that when 7 = 0, both walls 
can stably move with nonzero velocities and form stable 
stationary and oscillatory, quiescent and travelling bound 
states. The resulting bifurcation diagram will then be 
interpreted using two integrals of motion of eq.(^, and a 
relation of our model to biaxial ferromagnets established. 
Turning to the dissipative case, we will give an analytical 
proof of the stability of the Neel wall for all h and 7 and 
describe stable bound states formed by the damped walls. 
□ 

Both static kinks, (||) and (Q), belong to a broader 
class of uniformly moving solutions of the form Tp{x — 
Vt). We found these by solving equation ^tpxx — iVipx — 
IV'PV' + V' = hi})* [Q. Fig.|l|(a) shows the momentum 
of the travelling wall as a function of its velocity. The 
momentum P = Im J Tp^ip*dx is one of the two conserved 
quantities of eq.(^) with 7 = 0, and hence is a natural 
choice for the integral characteristic of its solutions. The 
second integral is energy, and it will also be used below: 
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dx. (5) 



(Here = 1 + h.) The stability of the travelling walls 
was examined |18 by computing eigenvalues A of 



(6) 



where the column ip — (u,v)'^, the operator TC is given 
by 



W = -7:91 
2 



37^2 +I^ + h 2711 -Vd.x+j 



and J is the antisymmetric matrix with J21 = —0^12 = 1- 
Eq. is obtained by linearising eq. (|^) about ip = TZ + iX 
in the co-moving frame, and letting Si}; = {u + iv)e^^ . 

In a striking contrast to the diffusive and relativistic 
dynamics |jll|, our numerical analysis of eq.(|^) 

reveals that not only the stationary Neel wall, but the 
entire branch of travelling kinks in Fig.0(a) is stable. This 
multistability admits a simple explanation in terms of 
the energy and momentum, though. The energy of the 
stationary Neel wall, Ejsj = |(1 + h)^^'^, is greater than 
that of the stationary Bloch wall, Eb = 4-\/7i— f ^'^^^i and 
so one might expect V'a' to decay into ipB plus radiation 
waves — as in the relativistic case |Q. However, unlike 
their relativistic counterparts, our Bloch and Neel walls 
have unequal momenta, with Pb > Pn (see Fig.|l|(a)) — 
and this makes the i})n ^ tpB decay impossible. 

Next, our simulations of the time-dependent eq.(|^) 
show that two stationary Bloch walls with opposite chi- 
ralities (i.e. opposite signs in (^) can attract and form a 
motionless breather-like bound state (Fig.|l|(b)). An at- 
traction of stationary Bloch and Neel walls results in a 
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FIG. 1: (a) The momentum of the travelling Bloch and Neel 
walls (thick) and their nonoscillatory bubble-like complex (thin 
line). For \V\ close to c, the wall attaches a small-amplitude 
bubble on each flank; this accounts for the turn of the thick 
curve near |V| = c. The dotted segments of the continuous 
branches indicate unstable solutions, (b) The formation of an 
oscillatory breather-like complex of two walls. (Only the real part 
of ip is shown for visual clarity.) In (a), h = y^; in (b), h = 0.1. 



moving breather. There also exist nonoscillatory, bubble- 
like, bound states of i})m and i/'s- Note that unlike their 
parent walls, all of these complexes approach the same 
background as a; ^ 00 and a:: ^ — 00. The bubble- like 
solitons admit most transparent physical interpretation: 
they describe "islands" of one stable phase in the sea of 
the other one, e.g. patches oscillating 180° out of phase 
with the rest of the vibrated water channel or chain of 
pendula. Below we focus on the bubbles and relegate the 
breathers to a separate publication. 

For each h there is a one-parameter family of motion- 
less bubbles, the parameter being the separation distance 
z between the two walls. (Accordingly, there are two zero 
eigenvalues in the spectrum of the operator ^ associ- 
ated with each bubble, one translational and the other 
one corresponding to variations in z). There is a partic- 
ular separation 2 = C foi' which the bubble is symmetric: 
ip^{—x) = —ipi^^x). The symmetric bubble turns out to 
have the largest momentum over bubbles with various z, 
while C is the smallest possible separation: z > C- More 
importantly, it is the only stable bubble. All nonsymmet- 
ric bubbles (z > C) were found to have a pair of nonzero 
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real eigenvalues ±A in their spectrum. As z ^ the 
pair converges at the origin and so the symmetric bubble 
has four zero eigenvalues. 

Let ±A be a pair of eigenvalues diverging from zero as 
z grows from ( and assume that A = e^/^Ai + . . . for small 
e — z ~ Then the associated eigenfunction ip expands 
as If = ipo + e^^^ifi + . . . , where (po is a linear combination 
of the two zero modes: tpo = (Cidxtp + C2dzip)\z=c, with 
■0 = (7^,I)^. (The other two zero eigenvalues have only 
generalised eigenvectors associated with them.) Substi- 
tuting this into (||) and using Tiz = Ti-c + eTii + . . . , the 
order e^/^ yields Ti(<-pi = XiJ<fo- A bounded <fi{x) exists 
only if XiJipo is orthogonal both to dz'>jj and dx4'- This or- 
thogonality condition amounts to {dP/dz)\z=(;; = 0, and 
the latter relation explains why the momentum has to 
reach its maximum at the value of z for which a pair of 
real eigenvalues converges at the origin. 

The other implication of the relation dP/ dz — is that 
it allows the symmetric bubble to be continued to 7^ 
[^9| . The resulting branch of moving bubbles is shown in 
Fig.|l|(a). As \V\ ^ c = {1 + 2h + ^Ah{l + h))^^^ , which 
is the minimum phase velocity of linear waves, the bub- 
ble degenerates into the flat background, whereas when 
y, P ^ 0, it transforms into a pair of Neel walls with the 
separation z — > 00. The entire branch of moving bubbles 
is stable, with the exception of a small region between 
V — and the point of the maximum \P\ inside which 
a real pair ±A occurs (FigJ^(a)). The change of stability 
at points where dP/dV = 0, is explained in jiot . 

The diagram Fig.|l|(a) can further justify referring to 
the kinks (||) and (jj) as Neel and Bloch walls. In fact sta- 
tionary nonchiral interfaces called Neel walls are known 
in uniaxial ferromagnets, where they coexist with chi- 
ral (Bloch) walls. When the axial symmetry is broken, 
the two types of walls can move; this occurs in partic- 
ular in easy-axis ferromagnets with the second, weaker, 
anisotropy axis (/3, e < and iJ = in eq.(||) below) 
p6[ . The P{V) curve for the easy-axis walls is qualita- 
tively similar to our Fig.|l|(a): The Neel wall's momentum 
Pn = while Pb 7^ 0; as the velocity V grows, the two 
branches are drawn closer together and finally merge. 
The limit velocity V — w is known as the Walker's ve- 
locity ll^. The E{V) curves are also similar. 

This analogy suggests that there could be a link be- 
tween the time-dependent NLS (||) and Heisenberg ferro- 
magnets and indeed, there is one. Consider a quasi-one- 
dimensional ferromagnet with a weakly anisotropic easy 
plane {Mx, My), in the external stationary magnetic field 
along Mz. The magnetisation vector M = {Mx, My, Mz) 
lies on the sphere, — Afp . and satisfies the (damped) 



Landau-Lifshitz equation |16 



-M^ =Mx [ W dC~ XMxMr, (7) 
oM J 

W = I (a^ M)2 + + '-^Ml - HMz + Wo, (8) 



2^, 



with /3 > 0. If the anisotropy parameter e is small and 
the field H is close to /SMq: H = (iM^ - eq, the vector M 
will stay near the northern pole of the sphere. Choosing 
s = qMo - PM^/2 > for e > and s < for e < 0, 
we define Mx + iMy — (2es//3)-^/^V-'*- Assuming that 
the relaxation constant A is ©(e^/^) or smaller, and that 
M depends only on "slow" variables x — (es/2aMQ )^/^^ 
and t — {2e^j,os/hMo)T, eq.(0)-(^ reduces to eq.(||) with 
h = (3M^/{2s) and 7 = 0. Note that the resulting NLS 
is undamped — although the original Landau-Lifshitz 
equation did include a small damping term. The effect 
of damping will become noticeable only on time scales 
longer than e~^; these are not captured by eq.(^. Also 
note that despite the analogy between the easy-axis and 
easy-plane ferromagnets, there are important physical 
differences. In particular in the easy-axis case the walls 
interpolate between M/Mq = i(0, 0, 1) while in our case 
they separate domains with M ~ (0,±|e|^/^, 1). □ 

Proceeding to the damped situation, 7 7^ 0, our first 
goal is to demonstrate the stability of the Neel wall. We 
let ip(x, t) — iPn{x) + S'ip{x, t), where 



5i}{x,t) = [u{X) + iv{X)]e^>'' 



T)T-ie. 



(9) 



X ~ Ax, T — A^t, V ~ A and ji is complex. Linearis- 
ing eq.(^ in small Sip we obtain an eigenvalue problem 



(Lo + e)w — {^^ - r)M, Liu = —{p. + T)v, 



(10) 



where e = 2 — 2/^4^ and Lq and Li are the Schrodinger 
operators with familiar spectral properties: 



1«2 



di^ - secliX, Li = Lq + 2 tanh^AT. 

Introducing = — and w = v^^{ijl + T)v ||], we 
eliminate F from the eigenvalue problem (|o|): 

{Lq + e)w — vu, Liu = —vw. (11) 

Now we will show that < for all < e < 2, so that 
^2 ^ p2 g^j^j-j perturbations decay to zero as i — > 00. 

The operator Li has a zero eigenvalue, with the asso- 
ciated eigenfunction yo{X) = sech^A, and no negative 
eigenvalues. Consequently, on the subspace TZ defined by 



u{X)y^{X)dX = 0, 



(12) 



there exists an inverse operator L^^ and so ( pH ) becomes 
(Lo + e)u — —v^L^^u, with Lo + e symmetric and a 
positive operator. The smallest eigenvalue —Vq is given 
by the minimum of the Rayleigh quotient: 



f u{Lq + e)udX 

mm - — . 

uen j uL^ '-udX 



(13) 



To prove that — t'g > it is sufficient to show that the 
minimum of the quadratic form J u{Lq + e)udX is pos- 
itive on TZ [EtI. Assuming that u{X) are normalised by 
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J v?dX = 1, the minimum is attained on the solution 
u{X) to the nonhomogeneous boundary- value problem 



{L^ + e)u{X)=riu{X)+ayo{,X), 



(14) 



where "q and a are the Lagrange multipliers. The mini- 
mum equals 7] — provided rj and a are chosen so that the 
u{X) satisfies eq.(|l2|) and the normalization constraint. 

The operator Lq has a single discrete eigenvalue Eq = 
— ^ with the eigenfunction zq[X) = (l/'\/2)sechX, and 
the continuous spectrum of eigenvalues E{k) = fc^, with 



Zk{X) 



ik + tanhX 
ik-1 



— ikx 



-oo < k < oo. 



(15) 



Expanding j/o and u over the complete set {zo] Zk} gives 



yoiX) = YoZoiX) + J Y{k)zk{X)dk, 
u{X) = UqZq{X) + I U{k)zk{X)dk. 



Substituting into (y_4|) and using the orthogonality of the 
functions in the set produces U{k) — alk"^ + e — rj)~^Y [k) 
and C/q = a{EQ + e — vi)^^Yt^. Using these in (f2|) gives 



^0 



En 



e-T] 



\Yik)\' 



-dk = 0. 



(16) 



The minimum of the quadratic form J u{Lq + e)udX 
is given by the smallest root ij* of the function (|l6| ) . The 
function is increasing for — oo < rj < e, apart from 

the point rj = Eq + c where it drops from +oo to — oo. As 
77 —00, Qeiv) ~^ +0; as 77 ^ e, 2e{v) tends to a finite 
value. (This follows from the fact that 



y(fc) 



ik 



TTk/2 



l + ik sinh(7rfc/2) ' 



hence the integral in (|l^) converges for all -q < e.) Con- 
sequently, there is only one root ij* and its sign is oppo- 
site to the sign of 0e(O). Since dQf{r])/de < 0, we have 
0e(O) < 00 (0) while the value 0o(O) can be calculated as 

0o(O) = J + ^^^fc = / yoLo'yodX. (17) 

Noticing that L^^yolX) = — 1 -I- ctanhX, with c an ar- 
bitrary constant, eq.(|T7|) yields 0o(O) = —2 and hence ij* 
cannot be negative for any e. Thus — i^^ > and the Neel 
wall is stable for all h and 7 (with h> j > 0). 

Are there any other attractors for nonzero 7? When 
7 7^ 0, the momentum is, in general, changing with time: 
P = — 27P, and therefore a uniformly moving soliton 
has to satisfy P = 0. Since both curves in Fig.|l|(a) cross 
the P = axis only at the stationary Neel walls, we 
conclude that no other solutions persist for small nonzero 
7. This does not, however, exclude the existence of new 



solutions for larger 7. Our numerical analysis has, in 
fact, revealed a window of h values, 7 < /i^ 0.35 -I- O.87 
(with 0.1 < 7 < 0.85) where two Neel walls attract and 
form a stable stationary bubble. □ 

In conclusion, the remarkable stability of the damped- 
driven kinks and their bound states is in sharp contrast 
with stability properties of the bright solitons. The sta- 
ble coexistence of two types of domain walls and their 
complexes in the undamped case is also worth emphasis- 
ing. This multistability is not observed in the paramet- 
rically driven Klein-Gordon and Ginzburg-Landau equa- 
tions and is due to the availability of the momentum in- 
tegral which takes different values on different solutions. 
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